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Abstract 

^  Both  one-dimensional  and  three-dimensional  models  are  presented 
for  the  collision  of  a  vibrationally-excited  molecule  with  a  semi¬ 
conductor  surface,  where  the  transfer  of  vibrational  energy  leads  to 
the  formation  of  electron-hole  pairs.  The  transition  probability  P  is 
calculated  as  a  function  of  the  molecule-surface  distance  for  real  two 
systems,  HC1  +  InSb  and  HC1  +  PbSe,  as  well  as  for  some  model  systems 
with  different  values  of  parameters.  While  P  generally  increases  as 
the  distance  decreases,  there  are  some  minima  at  intermediate  distances. 

The  overall  probability  is  obtained  as  an  approximate  integral  of  P  over  the 
distance,  and  for  thermal  collisions  values  of  a  few  percent  are  obtained. 
Such  values  are  high  enough  for  an  experimental  observation  of  electrical 
conductivity  due  to  electron-hole  pair  formation. 
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I.  Introduction 

A  vibrationally-excited  molecule  near  a  solid  surface  can  transfer  vibrational 
energy  to  either  phonons  or  electrons  (or  both)  of  the  solid.  The  situation  involving 

1  O 

excitation  of  phonons  has  been  actively  studied,  and  theoretical  investigations 
of  energy  transfer  to  the  electrons  of  a  solid  have  been  carried  out  for  the  case 
of  a  metallic  surface.  In  the  present  study  we  consider  the  transfer  of  vibrational 
energy  to  the  electrons  of  a  semiconductor  surface.  This  process  differs  from  the 
case  of  a  metallic  surface  due  to  the  energy  gap  between  the  valence  and  conduction 
bands.  Here  the  transfer  of  vibrational  energy  can  excite  electrons  from  the  valence 
to  the  conduction  band  and  hence  result  in  the  formation  of  electron-hole  pairs  which 
can  be  detected  by  measuring  the  semiconductor  electrical  conductivity.  *  In 
Section  II  we  present  a  one-dimensional  model  which  depends  on  a  transition  matrix 
element.  This  element  is  evaluated  in  Section  III,  and  the  results  are  used  to 
determine  the  transition  probability  in  Section  IV.  The  polarization  of  the  solid 
state  is  considered  in  Section  V,  and  an  extension  to  three  dimensions  is  discussed 
in  Section  VI.  The  results  of  numerical  calculations,  with  specific  applications  to 
the  HC1  +  InSb  and  HC1  +  pbSe  systems,  are  presented  in  Section  VII,  and  Section  VIII 
is  the  Summary  and  Concluding  Remarks. 
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II.  The  Model 

Let  us  consider  a  diatomic  molecule  whose  center  of  mass  is  located  at  some 
distance  R  from  a  semiconductor  surface.  Within  the  framework  of  a  one-dimensional 
model,  the  molecule  is  assumed  to  be  oriented  perpendicular  to  the  surface  on  the 
continuation  of  a  line  formed  by  a  linear  chain  of  lattice  atoms  (Fig.  1),  where  the 
interaction  of  the  molecule  with  the  solid  is  restricted  to  just  the  atoms  on  the 
chain.  We  assume  the  velocity  of  the  molecule  to  be  sufficiently  low  such  that  the 
influence  of  its  translational  motion  can  be  neglected.  While  the  vibrational 
motion  of  the  molecule  affects  the  motion  of  both  the  nuclei  and  electrons  of  the 
lattice,  we  shall  consider  its  effect  only  on  the  electrons. 

The  transfer  of  molecular  vibrational  energy  to  the  chain  electrons  leads  to  the 
excitation  of  electrons  from  the  valence  band  to  the  conduction  band  and  therefore 
to  the  formation  of  electron-hole  pairs.  Such  excitation  is  possible  only  when  the 
vibrational  energy  level  spacing  Ey  a'tfw  is  larger  than  the  semiconductor  band  gap 
energy  (Fig.  2), 


Ey  =irW>  Eg. 


(1) 


For  simplicity  we  take  the  molecule  to  be  in  its  first-excited  vibrational  state 
and  assume  the  influence  of  the  solid  on  its  vibrational  motion  to  be  negligible 
so  that  Ey  is  constant.  We  further  restrict  ourselves  to  cases  where  Ey  is  less 
than  2Eg  so  that  only  a  one-electron  transition  is  possible. 

Considering  the  coupling  between  the  vibrational  motion  of  the  molecule  and  the 
solid  electrons  as  a  weak  perturbation,  we  can  use  the  golden  rule  to  determine  the 
transition  probability,  which  can  be  written  for  the  case  of  a  one-electron  tran¬ 
sition  between  two  conti nua  as 

eT 


r  v 
,  —  I 

~  *tT  J  jnin 


^min  dEV  VEV}  nC^EC^  W01^EV,EC^ 


(2) 
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P  is  the  probability  of  transition  per  second,  Ey  and  Eg  are  one-electron  energies 
of  the  valence  and  conduction  bands,  respectively,  riy  is  the  density  of  electron 
states  in  the  valence  band,  nc  is  the  density  of  empty  levels  in  the  conduction 
band,  and  Wg1  is  the  matrix  element  for  the  transition  from  the  excited  vibrational 
state  to  the  ground  state  with  a  simultaneous  transition  from  the  valence  band  to 
the  conduction  band.  Due  to  energy  conversation,  for  this  transition  the  Eg  level 
is  greater  than  Ey  by  the  vibrational  quantum  Ey, 

Ec  *  Ev  +  Ev  •  (3) 

The  upper  limit  on  the  integral  is  the  maximum  electronic  level  of  the  valence  band, 
and  the  lower  limit  is  given  by 

Eyin  *  Ejax  -  (Ev  -  Eg)  .  (4) 

The  matrix  element  WQ1  is  expressed  as  an  integral  over  the  electronic  coordinate  • 

z  and  the  vibrational  coordinate  Z  (Fig.  2), 

W01(Ey,Ey+Eg)=JdzdZ  4^(z)  x^Z^z.Z^Ie  (zJxgtZ),  (5)  • 

where  4/^  and  4/^  are  the  electronic  warefunctions  in  the  valence  and  conduction 

bands,  respectively,  Xq  and  are  the  ground  and  first-excited  state  vibrational  • 

A 

warefunctions,  respectively,  and  W  is  the  interaction  between  the  electron  and  the 
molecule.  At  close  distances,  where  the  overlap  between  the  molecular  orbitals  and 

A 

the  electronic  states  of  the  solid  cannot  be  neglected,  W  is  complicated.  At  far  • 

enough  distances  where  the  overlap  can  be  neglected,  W  becomes  electrostatic.  If 

the  size  of  the  molecule  is  smaller  than  the  distance  to  the  solid-state  electron, 

then  the  field  of  the  vibrating  molecule  can  be  expressed  as  a  sum  of  dipole  and  • 

higher-order  multipole  potentials.  We  shall  consider  only  the  point  dipole  potential, 
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as  has  been  done  for  the  case  of  a  metal  surface.®  If  the  polarization  of  the  solid 
is  neglected,  i.e.,  setting  the  dielectric  constant  equal  to  one,  then  the  interaction 
becomes  (in  atomic  units). 


I 


d 


i--eUL  ,  (6) 

(R-z)2 

where  p(Z)  is  the  dipole  moment.  Substituting  Eq.  (6)  into  (5),  we  obtain  a  separa¬ 
tion  of  variables. 


”oi(EVEVtEv>  "  P01  /_”  dz  4y)(z)  ^2 

where  L  is  the  length  of  the  chain  and  Pq1 
1  -*■  0  transition. 


4CL  (z)  ,  (7) 

CV  cv 

is  the  dynamic  dipole  moment  for  the 
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Pqi  =  f  Xo(Z)p(Z)Xl(Z)  .  (8) 

The  separation  of  variables  is  also  achievable  when  other  multipole  moments  are 

8 

included  in  W.  The  dynamic  dipole  moment  is  known  for  many  diatomic  molecules  from 
experiments  and/or  ab  initio  calculations.  The  determination  of  is  thus  reduced 
to  the  evaluation  of  a  one-dimensional  integral  over  the  solid-state  electronic  coor 
dinate.  In  the  model  considered  here,  only  polar  molecules  can  transfer  vibrational 
energy  to  the  solid.  However,  at  close  distances,  particularly  for  an  adsorbed 
molecule,  this  transfer  is  possible  also  for  non-polar  molecules  due  to  the  overlap 

i 

between  the  molecular  orbitals  and  the  electronic  states  of  the  solid. 


) 


» 
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III.  Transition  Matrix  Element 

In  order  to  calculate  the  transition  matrix  element  displayed  in  Eq.  (7),  one 

must  know  the  solid-state  electronic  wavefunctions  and  Forms  for  these 

11  12 

are  available  in  terms  of  the  wavenumber  k,  but  we  can  reexpress  these  in 
unitless  and  more  symmetrical  forms  in  terms  of  the  parameter  x  determined  from  the 
equations 

k  =  f(l+x).  -K<Ul  (9) 

x  =  2n  .  i  =  2a  n.  i,  (10) 

where  g  =  Nk/n  =  2ir/a  is  the  reciprocal  lattice  constant,  a  is  . the  lattice  constant, 
n  is  an  integer  (0ini.N),  and  N  =  L/a  is  the  number  of  unit  cells.  The  wavefunctions 
and  also  the  energies  then  can  be  written  as 


a.  ■  i  pFJ  4V)(u)  +  Jffi  *<U>(u)] 

(ID 

*'C)<u)  *  -J-pFJ  *<c)(u)  +  /r?  ♦<c)t«)l 
x  ^  L  x  x  J 

(12) 

« 

Si  sin[ir(l+x)(u-ud)  +  ©x] ,  u<ud 

(13) 

*<V)(x)  ■  ■ 

A 

4 

Si  sin0x  expE-Tr/8r-(l+x)2(u-ud)  ]  ,  u>urf 

(14) 

Ey(x)  =  | 

[l+x2-  lU c2-fx2  J 

(15) 

4 

-Si  C0S|7t(1+X  )  (u-ud)  +0X  ,  U>Ud 

L  J 

(16) 

*<C)(u)  - 

i 

-si  cosexexp^-Tr/^r-(l+x)^(u-ud)j  ,  U£Ud 

(17) 

■ 

(17) 
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EC(X)  =  f-  [l+x2+2/4^+x2  , 


(18) 
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where  u  =  z/a,  =  d/a,  d  is  a  distance  on  the  order  of  an  atomic  radius  (an  actual 

2  2 

value  is  set  in  Section  VII),  x,  =  E^/g  ,  r  =  W/ g  ,  W  is  the  sum  of  the  work  function 
and  Fermi  energy,  y  =  t4e2+x*\  and 
tane  =  (l+x)/v4r-(l+x)2. 


(19) 


(V)  (Cl 

The  wavefunctions  45 v  '  and  <f>'  '  are  normalized  to  the  unit  length. 


Substituting  the  above  forms  of  the  wavefunctions  into  the  integral  of  Eq.  (7) 


and  extending  the  lower  limit  of  integration  to  we  obtain  for  arbitrary  4* 
and  : 


(V) 


Pm 

Wqi  =  —  Q(x  ,x) 


Q(x,x)  =  ~ if^y-xKy+x)  S  5  +  ✓(y-x)(y-x)  S  r 

2/yy  *-  x’x  x’_> 


+  /(y+x)(y+x)  S  ~  +  /(y+xKy'x)  S 

“  A  }  A  - 


-■*] 


(20) 


(21) 


fl 


S  -  =  S^1-  +  S(I*> 
x,x  x,x  x,x 


s[l)~  =  -2sinevcos0?  f  dt  exp[-rr(t4r-(l+x)2 

X,X  X  X  Jq 

+  v^r-(l+x)2)t] 


(p-t) 


2  ’ 


(22) 


(23) 


t  =  u-u. 


^  -  *r « °(t) 


(24) 


Q(t)  =  sin[ir(l+x)t  -  9  ]  COS  [it  ( l+x)t  —  0]  - ; 

x  x  (p+t>‘ 


(25) 


t  =  ud-u 
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where  p  =  ur  "  u(j  an^  uR  =  R/a. 

The  integral  in  Eq.  (23),  which  describes  the  interaction  with  the  decaying  por¬ 
tions  of  the  wavefunctions  outside  the  solid— Eqs.  (14)  and  (17),  goes  to  infinity 
at  t  =  p.  However,  for  t  close  to  p  the  point  dipole  approximation  is  not  valid 
and  the  finite  size  of  the  molecule  must  be  taken  into  account.  Fortunately,  the 
integral  in  Eq.  (23)  is  usually  smaller  than  the  integral  in  Eq.  (25)  so  that  the 
former  can  be  treated  approximately.  The  finite  size  of  the  dipole  can  be  accounted 
for  by  replacing  the  term  (p-t)  by  (p-t)  +  pM»  where  pm  is  a  distance  on  the  order 
of  the  relative  (per  lattice  constant  a)  size  of  the  dipole  (molecule)  This  proce¬ 
dure  is  valid  only  for  p  >  p^.  Usually  the  exponential  function  in  E  (23)  de¬ 
creases  quickly  and  becomes  small  when  t  >  p^  Accounting  for  this  I-  vior  and 

o 

taking  p  to  be  greater  than  pM>  we  may  neglect  the  t-dependence  of  the  term  (p-t)  + 

2 

PM  and  express  Eq.  (23)  roughly  as 


2sir.oxcose- 

ir(v&r-(l+x)^  +  \^5r-(l+x)2)(p2+pjj) 


(26) 


The  integral  of  Eq.  (24)  which  describes  the  interaction  with  the  bulk  portions 
of  the  wavefunctions  is  expressed  in  terms  of  sine  and  cosine  integrals  as 


S(II)  =  -T(  +  )  -  T(-} 

x,x  x,x  x,x 

Txt;  *  >»<  W  +  ^ 

tl  ■  ♦  *<*-*> 


(27) 

(28) 

(29) 


U 


(  +  ) 

x  ,x 


=  COS(ira+Q  +0' )ci  (to)  +  sin(7ra+9  +P;)si  (ttcx) 

A  A  A  A 
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(30) 
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=  cos(ire49x-e~)c  1(7:8)  +  si  n(  Trg+e  x~e^)si  (ne) 


(31) 


where  a.  =  (2+x+x)0  ?n-  p  =  (x-x)p. 
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IV.  Transition  Probability 

The  transition  probability  can  be  calculated  according  to  Eq.  (2)  by  integra¬ 
ting  over  the  electron  energy.  However,  it  is  more  convenient  to  perform  the 
integration  in  terms  of  the  variable  x  which  determines  the  wavenumber  as  in  Eq.  (9). 
The  integration  must  be  performed  from  -Xq  to  +Xq,  where  Xq  corresponds  to  the 
energy  Eyin,  also  taking  into  account  that  for  each  x  there  are  two  transitions: 
x  x  and  x  -*■  -x  with  transition  energy  equal  to  the  vibrational  energy  Ey  (Fig.  3). 
Due  to  the  symmetry  of  signs  in  Eq.  (21)  --  Q(x,x)  =  Q(x,-x)  =  Q(-x,x)  =  Q(-x,-x), 
we  can  change  the  lower  integration  limit  to  0  simply  by  multiplying  the  integral  in 
Eq.  (2)  by  4.  Substituting  Eq.  (21)  into  (2),  we  then  have 


where  Xy  corresponds  to  the  energy  Ey  of  the  valence  band  and  x^.  to  the  energy 
E^.  =  Ey  +  Ev  of  the  conduction  band,  and  tq  =  =2.419x10  ^sec  is  a  time 

atomic  unit  which  has  been  inserted  so  that  the  quantities  in  the  above  equation 
can  be  expressed  in  atomic  units,  giving  P  in  sec  . 

Each  level  of  the  valence  band  is  occupied  by  two  electrons  of  opposite  spin, 
so  that  riy  is  obtained  by  doubling  the  level  density,  and  Oq  is  simply  equal  to 
the  level  density  due  to  conservation  of  spin  in  the  band-to-band  transition: 


y(xy)  -  2 


dn^ 

dT 


dE, 


dnr  .  dxr 
r'C^xC^  =  dE^  =  2a  I  dE^ 


(33) 


(34) 


where  n  and  nr  are  the  numbers  of  levels  with  the  wavenumber  valuer  §  U+Xv)- 
Eq.  (10)  has  been  used  to  obtain  dn/dx  =  L/2a  =  l/2a,  where  the  same  normalization 
is  used  for  the  level  densities  as  for  the  wavefunctions  in  Eqs.  (13),  (14),  (16) 
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a 


and  (17)  which  were  normalized  to  the  unit  length. 

Combining  the  derivative  of  Eq.  (18)  with  Eqs.  ( 32 )- ( 34) ,  one  can  write 


P  = 


4A  /  P01 


f(uD) 


a  r'“R' 

X, 


f  ^XV  ^ 


/4^+xj 


xc(/4?2+x2-l) 


where 


=  ^  ^[l-2?(2£-l)]  - 

=  |[2+8E?-2y+x2]  -  2/+4c2+8^-2y+x2j*  , 


(35) 

(36) 


(37) 

(38) 


E  =  Ey/Eg  >  1,  and  ?  and  y  are  defined  just  after  Eq.  (18).  At  the  point  xy  =  xQ, 
xc  equals  zero  but  the  integral  is  finite, because  in  the  vicinity  of  Xy  =  xQ  the 
variable  x^  is  proportional  to  VxQ  -  Xy  The  integral  must  be  evaluated  numerically. 

We  should  mention  that  the  transition  probability  can  also  be  expressed  directly 
as  an  integral  over  energy  as  follows: 


(39) 


( 


1 


where  E  =  8Ey/g2,  Ej  =  1-4?(2E-1),  E2  =  1-4?,  and 

yy  =  |  2v4?^+u  -  1  -  u  |  *  (40) 

yc  =  |  2i4?2+8C?+u  -  1  -  u  -  8E?  |*  . 


(41) 


V.  Polarization  of  the  Solid  State 

The  expression  for  the  transition  matrix  element  in  Eq.  (20)  was  obtained  by 
the  assumption  that  the  electrical  field  [Eq.  (6)]  in  the  bulk  is  the  same  as  in 
the  vacuum.  In  other  words,  the  dielectric  constant  e,  which  is  often  large  for 
semiconductors,  was  set  equal  to  one.  The  dielectric  constant  cannot  be  taken 
into  account  simply  by  introducing  it  into  the  interaction  in  Eq.  (6).  In  the 
vicinity  of  the  surface  the  dielectric  constant  is  a  function  of  the  distance 
from  the  surface,  changing  from  1  outside  the  solid  to  e  deep  inside  the  solid. 

It  is  just  this  region  at  the  surface  which  is  the  most  important  for  the  transition 
matrix  element  of  Eq.  (20).  Unfortunately,  the  introduction  of  any  reasonably 
smooth  function  e'(z)  into  the  integral  of  Eq.  (24)  precludes  an  analytical  solution. 
The  simplest  way  to  account  for  the  polarization  of  the  solid  is  to  make  e'(z)  a 
step  function, 

f  1,  z  ^  -  b 

e'(z)  •{  (42) 

,  z  <  -  b 

where  the  parameter  b  is  expected  to  be  on  the  order  of  the  diameter  of  the  Onsager 
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cavity,  which  can  be  assumed  to  be  the  diameter  of  an  atom.  Introducing  this 
into  Eq.  (24),  we  obtain 

=*7  f  dufi(u)  -  2fUd  Ufedufi(u)  ,  (43) 

Jud+ub  J° 

where  fi(u)  is  the  function  of  Eq.  (25)  and  =  b/a.  Evaluating  the  above  two 
integrals  in  the  same  way  as  for  Eq.  (24),  we  obtain 


c(II) 

°X,X 


(  +  ) 

x,x 


+  F 


-) 

,x 


» 


(44) 


where  the  first  two  terms  on  the  right-hand  side  are  given  by  Eqs.  (28)  and  (29). 
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The  last  two  terms  are: 


-p  sin[6x  +  9X  -  fr(ab-a)3  +  tt(2+x+x)Gx*~ 
pb 

i  sin[9x  -  6S  -  *(6b-S)]  +  •(*-«)«<'] 


GX^X  =  C0S^7ra  +  ex  +  ex}  +  sin(ra  +  6x  +  ex)  si  (TTOtfj) 

Gx"J  =  COs(itB  +  0X  -  8~)  ci(TTBb)  +  s i n ( ttB  +  9X  -  e~)  si(rrBb) 


where  =  (2  +  x  +  x)ub,  Bb  =  (x  -  x)ub  and  pb  =  uR  +  ub< 
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VI.  Three  Dimensions 

The  vibrationally-excited  molecule  can  obviously  form  electron-hole  pairs  not 
only  in  the  chain  on  the  same  line  with  the  molecule  (Fig.  1)  but  also  in  other 
chains  and  in  different  directions  of  the  electron  motion.  The  consideration  of 
this  three-dimensional  interaction  is  very  complicated,  and  in  order  to  simplify 
the  problem  we  shall  introduce  some  assumptions  which  enable  us  to  reduce  the 
three-dimensional  integration  to  an  integration  over  the  surface. 

The  solid  state  is  envisioned  as  formed  by  a  series  of  one-dimensional  chains 
which  are  parallel  to  each  other  and  perpendicular  to  the  surface.  The  interaction 
between  the  molecule  and  a  chain  which  is  not  on  the  same  line  as  the  molecule  is 
assumed  to  be  determined  by  the  expression  for  the  one-dimensional  case,  taking 
into  account  only  the  difference  in  the  distance  and  the  influence  of  the  dipole 
orientation  on  the  perpendicular  component  of  the  field  at  the  surface  point  A 
(Fig.  4).  Taking  the  lattice  to  be  formed  by  two  kinds  of  atoms  with  similar  sizes 
and  the  distribution  of  atoms  to  be  uniform,  we  can  then  express  the  transition 
probability  as 


P 


f(uR)  +  8ttI  du 


o 

"/ 


0.5 


urftuA) 


(49) 


where  ur  =  r/a,  u^  =  R^/a  =  +  u^,  and  is  the  distance  between  the  dipole  and 

the  point  A.  The  lower  limit  of  integration  is  taken  to  be  half  of  the  lattice 
constant,  considering  it  approximately  as  the  distance  from  the  point  0  to  the 
closest  neighboring  atom.  The  factor  (2uR  -  u^]?/4uj|  is  equal  to  the  square  of  the 
ratio  between  the  z-component  of  the  dipole  field  at  the  point  A  and  the  field  in 
the  one-dimensional  case  for  the  same  distance  from  the  dipole. 

The  above  approximations  better  describe  the  interaction  with  the  solid-state 
electrons  which  are  close  to  the  surface  than  with  the  electrons  deep  in  the  bulk. 
Fortunately,  when  the  dielectric  constant  is  large  the  electrons  deep  in  the  bulk 
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interact  relatively  weakly  with  the  molecule  and  hence  do  not  have  a  significant 
effect  on  the  transition  probability. 

< 
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VII.  Numerical  Results 

All  results  were  obtained  for  a  simple  lattice  of  two  kinds  of  atoms  (Figs.  1 
and  4),  where  the  one-dimensional  unit  cell  consists  of  two  atoms  only.  With  the 
assumption  that  the  atoms  do  not  differ  much  in  their  size,  the  parameter  d  [mentioned 
just  after  Eq.  (18)]  can  be  taken  as  equal  to  the  average  atomic  radius,  which  is 
approximately  a  quarter  of  the  lattice  constant,  i.e.,  a/4  (urf  =  0.25).  The  para¬ 
meter  b  for  the  dielectric  constant  in  Eq.  (42)  is  assumed  to  be  equal  to  the  atomic 
diameter,  a/2  (ub  =  0.5). 

The  transition  probability  given  by  Eq.  (35)  is  proportional  to  the  square  of 
the  ratio  Pg^/a,  where  Pg^  is  the  dynamic  dipole  moment  of  the  molecule  for  the 
1  ->  0  vibrational  transition.  The  dependence  on  the  molecule-surface  distance  as 
well  as  on  other  parameters  is  included  in  the  function  f(uR).  This  function  depends 
on  unitless  quantities  only,  such  as  the  relative  molecule-surface  dustance  uD 
[defined  just  after  Eq.  (25)],  the  relative  energy  gap  c  [after  Eq.  (18)],  the  sum  W 
of  the  work  function  and  Fermi  energy,  and  the  ratio  £  between  the  vibrational 
energy  and  the  energy  gap  [after  Eq.  (38)].  The  function  f(uR)  depends  very  weakly 
on  the  parameter  r ,  which  is  determined  by  W  [after  Eq.  (18)].  In  all  calculations 


we  take  r  to  be  0.7,  which  corresponds  approximately  to  W  =  5  eV  for  the  semicon¬ 
ductor  InSb. 

In  order  to  investigate  the  behavior  of  the  transition  probability  on  uR  as 

determined  by  f(uR),  we  have  carried  out  calculations  for  different  values  of  the 

parameters  £  and  £  and  the  dielectric  constant  e.  All  results,  except  for  the  one 

case  of  the  HC1  +  PbSe  system,  are  presented  for  the  lattice  constant  of  the  semi- 
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conductor  InSb,  a  =  6.48  A  (12.24  a.u.),  and  the  dynamic  dipole  moment  of  the 
HC1  molecule,  Pg^  =  0.0712  D  (0.0279  a.u.).^  The  results  for  other  values  of  Pg^ 
and  a  can  be  obtained  simply  from  the  equation 


P(ur;£,£,c) 


P(uR; £,£,e) , 


(50) 
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where  Pqj  and  a  are  expressed  in  the  units  D  and  A,  respectively. 

Results  for  the  dielectric  constant  e  =  1  are  presented  in  Fib.  5.  The  most 
interesting  feature  of  the  curves  is  that  the  transition  probability  P  does  not 
decrease  monotically  as  the  molecule-surface  distance  increases.  The  curves  (dashed 
lines)  for  the  one-dimensional  (ID)  calculation  exhibit  a  deep  minimum  which  is  shifted 
toward  the  surface  as  the  vibrational  energy  decreases,  i.e.,  as  £  decreases.  In  the 
three-dimensional  (3D)  calculation,  the  curves  (solid  lines)  are  smoother,  as  might 
be  expected  due  to  a  "smearing-out"  effect,  but  also  have  a  minimum  which  becomes  very 
shallow  for  large  vibrational  energies  (£  =  2.0).  At  large  distances  the  ID  tran¬ 
sition  probability  has  a  sharp  decrease  approximately  proportional  to  1/Rn  (n  =  2-4). 
The  decrease  is  slower  in  the  3D  case,  and  in  fact,  the  decrease  for  ^  =  1.1  does 
not  occur  until  uR  =  12.  [We  note  that  the  3D  results  obtained  by  the  integration 
over  the  surface  in  Eq.  (49)  are  very  approximate  for  c  =  1  when  electrons  deep 
in  the  bulk  of  the  semiconductor  play  a  significant  role.]  At  small  distances 
P  increases  quite  sharply  with  little  distinction  between  different  parameters  and 
models  (ID  and  3D).  Very  close  to  the  surface,  when  the  size  of  the  molecule  is 
comparable  to  its  distance  from  the  surface,  the  point  dipole  approximation  does  not 
work  well.  However,  for  a  relative  distance  uR  as  small  as  0.5  in  the  HC1  +  InSb 

o 

system,  for  example,  the  molecule  is  3.3  A  from  the  surface,  which  is  over  2.5  times 
as  large  as  the  size  of  the  molecule  and  hence  large  enough  for  the  point  dipole 
approximation  to  be  valid. 

Results  for  the  dielective  constant  e  =  15  are  presented  in  Fig.  6.  These 
curves  are  smoother  than  those  for  e  =  1,  and  for  small  vibrational  energy 
(E,  =  1.1)  they  decrease  monotonically  as  the  distance  increases.  Except  for  very 
small  distances,  the  transition  probability  is  smaller  and  decreases  more  steeply 
at  large  distances  for  e  =  15  than  for  e  =  1.  The  influence  of  the  dielectric 
constant  on  P  is  further  demonstrated  in  Fig.  7,  where  results  are  presented  for 
different  c  but  for  the  same  values  of  the  parameters  c  and  £;.  While  P  decreases 
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at  large  distances  as  e  increases,  the  opposite  behavior  occurs  in  some  intervals 
at  intermediate  distances. 

The  minimum  in  the  transition  probability  at  intermediate  molecule-surface 
distances  occurs  because  the  contributions  from  different  parts  of  the  chain  to 
the  integrals  of  Eqs.  (24)  and  (43)  have  different  signs.  This  leads  to  an 
interference  effect  which  is  less  important  for  large  e,  where  one  unit  cell  of 
the  chain,  namely  the  closest  to  the  surface,  makes  the  main  contribution  to  the 
integral . 

The  dependence  of  the  transition  probability  on  the  molecule-surface  distance 
for  two  real  systems,  HC1  +  InSb  and  HC1  +  PbSe,  is  presented  in  Fig.  8.  In 
contrast  to  the  results  discussed  earlier,  the  distance  here  is  not  in  relative 

O 

unitless  numbers  but  rather  in  A.  The  lattice  parameters  for  the  two  semiconduc- 
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tors  are  similar  --  a  =  6.48  A  for  InSb  and  a  =  6.124  A  for  PbSe.  The  main 
differences  are  in  the  energy  gap,  which  is  larger  for  PbSe  (0.26  eV)  than  for 
InSb  (0.165  eV),  and  in  the  dielectric  constant,  which  is  much  larger  for  PbSe  (280) 
than  for  InSb  (18). ^  At  large  distances  P  is  much  lower  for  PbSe  mainly  due  to 

o 

the  larger  value  of  e,  although  at  small  and  intermediate  distances  (4-21  A)  P 
is  higher  for  PbSe. 

It  is  possible  to  calculate  the  overall  transition  probability  corresponding 
to  each  curve  P(uR).  If  the  molecule  is  moving  with  velocity  v  perpendicular  to 
the  surface,  then  this  overall  probability  P^  is 

00 

PE=f  /  du  P(u),  (51) 

us 

where  u  =  R/a  and  u$  is  the  minimal  possible  distance  between  the  molecule  and 
surface.  It  is  reasonable  to  consider  this  to  be  the  distance  associated  with  an 
adsorbed  molecule.  However,  the  results  obtained  here  cannot  be  continued  to  such 
short  distances  due  to  the  approximations  of  the  model.  Nevertheless,  we  have 
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evaluated  the  above  integral  for  values  of  us  for  which  our  model  is  valid,  and  P^. 
is  given  for  three  different  values  of  u  in  Table  I.  The  dependence  of  P  on  the 

^  Li 

parameters  e,  x,  and  £  is  presented  in  Table  II  for  u$  =  0.5  and  0.9.  In  these 
calculations  the  3D  form  of  P(u)  was  used,  and  the  HC1  molecule  was  taken  to  be 
thermal  (300  K)  with  a  velocity  v  =  4. 0  x  10^  cm/sec.  We  see  that  P^,  is  on  the 
order  of  ten  percent  for  u  =0.5  and  is  on  the  order  of  one  percent  for  u  as  large 
as  0.9.  For  us=0.5  the  probability  does  not  depend  as  much  on  the  dielectric  constant, 
while  for  u$=  0.9  it  is  much  higher  for  e  =  1  than  for  e  =  15.  It  is  interesting  to 
note  that  as  the  vibrational  energy  parameter  £  increases,  P^,  decreases  in  spite  of 
the  increase  in  the  number  of  band  levels  involved  in  the  transition.  This  is  due 
to  the  fact  that  for  values  of  the  parameter  £  close  to  1,  the  transition  takes 
place  between  levels  whose  spacing  is  close  to  the  energy  gap,  so  that  the  product 
of  the  level  densities,  in  the  integral  of  Eq.  (2)  is  large  as  and  nc  are 

going  to  infinity  in  the  vicinity  of  the  gap. 
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VIII.  Summary  and  Concluding  Remarks 

For  a  vibrationally-excited  molecule  not  very  close  to  the  semiconductor 
surface,  the  influence  of  the  vibrational  motion  on  the  solid-state  electrons  can 
be  taken  as  an  electrostatic  interaction  with  a  point  dipole.  If  the  vibrational 
energy  level  spacing  is  larger  than  the  semiconductor  energy  gap,  this  interaction 
leads  to  the  excitation  of  electrons  from  the  valence  to  the  conduction  band  and 
hence  to  the  formation  of  electron-hole  pairs  which  can  be  detected  by  measuring 
the  semiconductor  electrical  conductivity. 

In  the  one-dimensional  model  the  transition  probability  (per  second)  is 
expressed  in  terms  of  integrals  of  analytical  functions  which  include  such  special 
functions  as  the  sine  and  cosine  integrals.  Results  from  this  model  show  the 
transition  probability  to  increase  sharply  at  short  distances.  At  intermediate 
distances  there  is  often  a  minimum  followed  by  a  maximum  and  then  a  further 
decrease  of  the  probability.  A  three-dimensional  calculation  is  performed 
approximately  by  the  integration  of  a  modified  one-dimensional  probability  over 
the  surface.  In  this  case  the  probability  is  higher  and  the  dependence  on  the 
molecule-surface  distance  is  smoother  than  in  the  one-dimensional  case.  The 
overall  probability  of  formation  of  an  electron-hole  pair,  obtained  by  summing 
from  the  distance  u$  to  infinity,  is  on  the  order  of  ten  percent  when  us  is  chosen 
to  be  0.5  of  the  lattice  constant  a. 

In  an  experimental  situation,  those  molecules  still  excited  at  u$  will 
continue  on  toward  the  surface  to  undergo  an  elastic  or  inelastic  encounter. 

Inelastic  processes  may  involve  either  phonon  or  electron  excitation  or  both. 

Because  the  events  occurring  after  the  molecule  reaches  u$  can  only  add  to  the 
probability  of  electron-hole  formation,  this  theory  provides  a  lower  limit  to 
the  efficiency  of  that  process.  Experimentally,^  efficiencies  in  the  range  0.1 

« 

to  10%  would  be  readily  observable.  In  summary,  this  model  calculation  points 
out  that  vibrational-to-electronic  energy  transfer  could  be  a  significant  energy 
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pathway  at  semiconductor  surfaces.  Possible  applications  include  selective  detection 
of  vibrationally  excited  molecules. 
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Table  I.  The  overall  transition  probability  (in  percent)  for  the  systems  HC1  +  InSb 
and  HC1  +  PbSe.  The  HC1  molecule  is  moving  perpendicular  toward  the  surface 
with  the  velocity  v  =  4.0  x  10^  cm/sec  from  infinity  to  the  point  u$  =  R/a . 


Table  II.  The  same  as  Table  I  but  for  HC1  +  semiconductor  with  different  values 
of  the  parameters  e,  c  and  £.  P^,  (0.5)  and  P^  (0.9)  are  given  in 
percent  for  u$  =  0.5  and  u$  =  0.9,  respectively. 
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Figure  Captions 

1.  Location  of  the  molecule  with  respect  to  the  solid-state  atoms  in  a  one¬ 
dimensional  chain;  a  --  lattice  constant,  R  --  molecule-surface  distance, 
d  --  border  distance  for  the  wavefunctions  of  Eqs.  (11)  and  (12). 

2.  Band  structure  of  the  semiconductor;  Z  --  vibrational  coordinate,  z  -- 
electronic  coordinate,  Ey3X  --  top  level  of  the  valence  band;  E^in  -- 
bottom  level  of  the  conduction  band;  Ey  --  vibrational  energy  level  spacing; 

Ey  &  E^  —  levels  involved  in  the  transition. 

3.  Band  energies  of  the  one-dimensional  chain  for  the  PbSe  semiconductor.  For 
the  valence  and  conduction  bands  Ey  and  Ec  see  Eqs.  (15)  and  (18),  and  for 
the  parameter  x  see  Eqs.  (9)  and  (10).  The  parameters  of  the  chain  are: 

a  =  6.124  A,  Eg  =  0.26  eV,  c  =  0.0324. 

4.  Three-dimensional  model;  0  --  center,  M  --  molecule,  R  --  molecule-surface 
distance,  r  --  distance  on  the  surface  from  0  to  atom  A,  R^  --  molecule-atom-A 
distance. 

5.  Dependence  of  the  transition  probability  P  (per  second)  on  the  relative 
molecule-surface  distance  uR  =  R/a.  The  value  of  the  dielectric  constant 

e  is  1.  The  numbers  on  the  curves  indicate  the  ratio  £  =  Ev/Eg  between  the 
vibrational  energy  level  spacing  and  the  energy  gap.  The  dashed  lines 
represent  the  results  of  the  ID  calculation,  and  the  solid  lines  represent 

p 

the  results  of  the  3D  calculation,  (a)  Relative  energy  gap  c  =  Eg/g  =  0.025 
and  (b)  t,  =  0.10.  [The  small  gaps  at  uR  =  5  for  the  3D  curves  are  because 
of  a  change  in  the  integration  stepsize  used  in  the  evaluation  of  the  integral 
in  Eq.  (49).] 


6.  Same  as  for  Fig.  5  but  for  e  =  15. 


7.  Dependence  of  the  transition  probability  P  (per  second)  on  the  relative  molecul 
surface  distance  uR,  where  c  =  0.05,  £  =  1.4,  and  the  numbers  on  the  curves 
indicate  the  dielectric  constant  e  (see  also  the  caption  to  Fig.  5). 

8.  Dependence  of  the  transition  probability  P  (per  second)  on  the  molecule-surface 
distance  R  for  the  systems  HC1  +  InSb  and  HCl  +  PbSe.  The  vibrational  energy 
level  spacing  is  0.358  eV.  The  parameters  for  HCl  +  InSb  are:  e  =  18, 

O 

a  =  6.48  A,  Eg  =  0.165  eV,  £  =  0.023,  £  =  2.17.  The  parameters  for 
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